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Abstract 

Using 586 pb^"*^ of e~^e~ collision data acquired at ^/s = 4.170 GeV with the CLEO-c detector 
at the Cornell Electron Storage Ring, we report the first observation of D^^ — t- D+e+e^ with 
a significance of 5.3(7. The ratio of branching fractions i3(£)*+ Dfe'^e^)/B{D*^ — ^ Dfj) is 
measured to be [0.72lQ;];3(stat) ± 0.10(syst)]%, which is consistent with theoretical expectations. 
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Dalitz decays in which a virtual photon is internally converted to an e~^e~ pair, 
have been observed in several vector-to-pseudoscalar decays of light mesons {e.g., u — )■ 
7r°e"'"e~, — )• 7r°e"'"e~, and — )■ 'r]e~^e~) However, such decays have not been reported in 
electromagnetic decays of mesons containing charm or bottom quarks. This Letter reports 
the first observation of such a decay, D*^ — t- D^e^e~ , and a measurement of its branching 
fraction. Only two decay modes of the D*~^ have been previously observed, the dominant 
£)*+ _^ mode and the isospin-violating D*"*" — )■ D^ir^ ^] decay. Their branching 

fractions have been determined by the PDG [2] from measurements of the ratio B{D*'^ — >■ 
D^'K^)/B{D*^ — > -Dj"7) and the assumption that they are the only D*^ decay modes. 

The expected D*^ Dalitz decay rate may be calculated by treating the photon from 
D*~^ — )■ as virtual and coupling it to an e~^e~ pair. The g^-derivative of the ratio 



can be written as 
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where q is the four-momentum of the virtual photon, rrix represents the mass of particle 



X, 



A = 



ml. 



— fnj^^, and fiq"^) is the transition form factor for D*'^ to . Motivated 
we use f{q'^)/f{0) = (1 — q'^/m'^)~^. Integrating Eq. ([2]), we 



by vector-meson dominance, 
predict Ree = 0.65%. 

We use ~ 5.6 x 10^ e+e~ — )■ DfO^^ events obtained from 586 pb'^ of e+e^ collision data 
with ^/s near 4.170 GeV acquired by the CLEO-c detector at the Cornell Electron Storage 
Ring (CESR). The CLEO-c detector is equipped with a CsI(Tl) calorimeter jEj to detect 
photons and determine their directions and energies, and two concentric cylindrical wire drift 
chambers [6| to track the trajectory of charged particles. The tracking chambers operate 
in an axial 1 T magnetic field to provide momentum measurements. The beam pipe and 
the drift chambers present under 2% of a radiation length of material, minimizing multiple 
scattering of charged particles and photon conversions. Charged hadron identification is 
achieved using energy loss {dE/dx) in the drift chambers and Cherenkov radiation in the 
RICH detector 0, [J. 

The default Kalman filter track reconstruction used to process CLEO data includes cor- 
rections for dE/dx and multiple scattering in the beam pipe and detector material, assuming 
each track has the mass of a pion, kaon, and proton. The tracks in this analysis are rather 
soft, with energies below 150 MeV, where dE/dx is very different from that of any of those 
three mass hypotheses. Therefore, to improve sensitivity, we reprocess events containing at 
least one exclusively reconstructed candidate, adding an mass hypothesis for each 
charged particle. Details of this analysis appear in Ref. |9|. CLEO has previously observed 



two other Dalitz decays, rj — )■ e~^e 7 l3] and 



7] — )■ e^e p 



111, in both of which the were 



substantially more energetic and did not need reprocessing. 

We reconstruct D*~^ — )■ Dfe~^e~ candidates using the nine distinct hadronic decay modes 
of the listed in Table [H Charge conjugate modes are also included. Candidates for Kg 



3 



and 7] are reconstructed through their decays to vr+vr" and 77, respectively. Measurement 
of Ree instead of the absolute branching fraction — )■ Z)^e+e~) bypasses any need 
for estimating the total number of D*~^ produced and minimizes systematic uncertainties 
stemming from reconstruction of the D^. 

We follow a blind analysis procedure to avoid bias. Selection criteria are optimized indi- 
vidually in each of the D'^ decay modes for maximum signal significance using Monte Carlo 
simulated samples of signal and background processes. The decay chain of signal events 
e~^e~ — )■ D*~^D~;D*~^ — )■ Dfe~^e~ are simulated with full angular correlations. Simulated 
samples of all e^e^ — )■ qq processes at 4.170 GeV where q = u,d, s, or c are used for back- 
ground. The reconstructed e"'"e~ tracks are required to pass within 5 cm of the interaction 
point in the direction parallel to the beam-axis and within 5 mm of the beam-axis in the 
transverse directions. The dE/dx of each candidate is required to be within So" of that 
expected for electrons. All charged pions and kaons in the decay chain are identified 
as such using a combination of dE/dx and RICH information as described in Ref. T2|. We 



require the reconstructed Df mass Mr,^ to be within a mode-dependent region around the 
known Dt mass l2| consistent with the resolution of the detector. We define the beam- 



constrained mass of the D*"*" by Mbc = y E'j^, — p|,, , where is the energy of the 
calculated from the beam energy and p^j is the three- momentum of the D*'^ inferred from 
its decay daughters' measured momenta. We select events with Mbc and 5M = Me>* — Md^ 
consistent with the known and masses [2|. 

A significant background to the observation of D*'^ D~^e^e~ arises from D*+ — t- -D^7 
events where the 7 converts into an e^e~ pair in the material of the beam-pipe or drift 
chambers. We reject much of this background using the following criteria for the tracks. 
We define the c/q of a track as the distance of closest approach of the track to the beam 
axis. Its sign depends on the charge of the track and whether the origin of the x — y plane 
falls within the circle of the track in that plane. We require the difference between the d^ 
of the e'^ and e~ tracks, Acio = d^ — d^ , to exceed —5 mm. Denoting each electron track's 
azimuthal angle measured at the point of closest approach to the beam axis by we also 
require A0o = 0o ~ 0o < 0.12. 

These selection criteria are applied on simulated samples of our signal to obtain the 
selection efficiencies for signal events e^g, where i stands for one of the nine decay modes of 
the used in this analysis. These criteria are applied to data in order to obtain the yields 
of events n'g. These numbers are presented in Tabled for each decay mode of the Z)+. 

Having established selection criteria using simulations, the background in the signal 
region for each mode i is estimated from the fit of an Mbc background function to the 
data in the Mbc sidebands for that mode. The shape of the Mbc function is fixed and is 
common to all modes. The function incorporates the kinematic limit and its parameters are 
determined from simulations. This shape is illustrated in Fig. 1(a). A similar estimate of 
obtained from the 5M distribution and the difference between this estimate and the Mbc 
estimate is taken as the systematic uncertainty inherent in the procedure. The estimated 
background for each mode is presented in Table |T] as 6gg. 

We calculate the signal significance, expressed in number of Gaussian standard deviations, 
from the Poisson probability for the estimated background to fiuctuate up to the observed 
signal yield or higher. The uncertainties in the background, both statistical and systematic, 
are modeled as Gaussian distributions. The combined signal significance for all modes is 
5.3(T. The most significant individual mode is K^K'tt^ at 5.0a. 

The distribution of the e~^e~ invariant mass Mgg for the 51 observed events in the signal 
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FIG. 1. Distributions of (a) Mbc and (b) 5M in data and simulated samples summed over all nine 
Df decay modes used in this analysis. In each figure, the points with error bars are data, and 
the unshaded histogram is the simulated D*g'^ — )• Dfe~^e~ signal. Background events in the upper 
shaded histogram (yellow online) are from simulated D*^ — t- Df^ decays. Background events 
in the lower shaded histogram (green online) are from simulated qq events that do not include 
D*'^ — )■ -0^7 or D*^ —7- Dfe'^e^. The curves are the fits of data to background shapes described 
in the text. The regions 2.100 to 2.124 GeV in (a), and 0.1298 to 0.1578 GeV in (b) are avoided 
in the shape fits to prevent contamination of the background estimates with signal events. 



region is compared to that expected in our simulations and presented in Fig. [2l The distri- 
butions are found to be in good agreement, with a Kolmogorov-Smirnov probability of 0.86 
for the events to share the same parent distribution. 

The criteria for selecting D*^ — )■ D^j events follow those for D*^ — t- Dfe^e^ as closely 
as possible. Instead of an e~^e~ pair, a photon candidate in the kinematically allowed energy 
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TABLE I. The yields n|e, estimated backgrounds 5gg, background-subtracted yields yjy, ratios of 
detection efficiencies and values R^^ of R^e for each D'^ decay mode i. The uncertainties given 
for ^* = e!y/egg are statistical only. The values of ngg and 6gg are summed over all modes, while the 
ratio of branching fractions i?gg is computed using Eq. ([3]) . 

i Ke ^ Ke 

K+K-TT+ 14 1.05t5]j^ ±0.79 9114 ±110 ± 201 4.65 ±0.12 0.66l[];^^ 

K+K-T:+n^ 6 1.70l[}:^^ ± 0.56 3592 ±118 ±72 4.80 ± 0.21 0.58t|j:^^ 

KlK+ 1 0.85i[]J^ ± 0.74 1902 ± 57 ± 45 4.31 ±0.13 0.03l[j:f^ 

K^gK'TT+TT+ 4 1.58t[}:|? ± 0-40 1570 ± 74 ± 13 5.38 ± 0.20 0.83l[]J|j 

7r+7r-7r+ 7 1.57t[}:^? ± 0.59 2745 ± 93 ± 52 4.62 ±0.10 0.9ltgjJ 

r/7r+ 4 1.40j:[}|^ ± 0.49 1037 ± 46 ± 37 3.87 ±0.10 0.97l[j:^^ 

r/p+ 7 2.62t[||| ± 0.23 3170 ± 161 ± 313 5.82 ± 0.24 0.80l[]j^ 

r/'7r+;??' ^ vr+vr-T? 4 0.00t[j:^g ± 0.00 691 ± 34 ± 40 3.96 ±0.12 2.30lJJ? 

??'7r+;7/' ^ /7 4 1.84^0:45 ± 0.25 1531 ± 80 ± 122 4.97 ±0.14 0.70i[];^^ 

All Modes 51 12. ei^};^^ ± 4.05 0-72^o:}3 



range is required and its electromagnetic shower in required to have a lateral spread narrower 
than 99% of true photons. The allowed window for 5M is broadened to account for worse 
resolution in calorimetric photon energy as compared to tracking-based e^e~ measurements. 
For similar reasons, the mass window restriction on Mbc is removed and the signal yield 
is extracted from a fit using background and signal shapes determined from simulations, as 
depicted in Fig. Ofor the K^K~i{^ mode. Application of these criteria to simulated samples 
of -D*'^ — 7- D^-^ gives us the efficiencies e!^. Application to data gives us the background- 
subtracted signal yields yl^ as listed in Table [B 

The ratio = e!y/e*g of efficiencies for D*+ D+7 to D*+ -> D+e+e~, given in Table [H 
cluster around 4.8. The K^K~ti^ mode stands out due to its large branching fraction, 
presence of two charged kaons and the absence of photons in the final state. No other 
mode or combination of modes is predicted to have more significance for an observation 
of D*'^ — )■ D~^e^e~ . For — i- K^K't^^ (all other modes combined), we expect a yield 
of ~15 (36) events, about 1 (13) of which are background. After unblinding the data as 
shown in Table [H we find yields quite close to those expected from our predicted i?eg 
and Monte Carlo simulations. Yields from individual non-K^ K~ t^^ modes range from 1 to 
7, consistent with expectations. 

For each of the decay modes i, we calculate the value of i?gg = (rigg — &ee)/(l/^/0 as 
listed in Table [B Note that the values of ^* are all equal (~ 4.8) to within ±20% because the 
dominant difference in efficiency is due to the photon versus e^e~ selection criteria, which 
are identical for all modes. The ^* vary somewhat with i due to the broadened 5M windows 
and signal shape Mbc fit (instead of a fixed window) for the radiative Dl^ modes relative 
to those of the Dalitz decays. Hence most systematic uncertainties in ^* due to selection 
cancel; e.g. those due to track-finding, particle identification, and selection of photons, vr", 
decays, etc. In order to avoid a bias due to Poisson fiuctuations when n^g is small and to 
preserve the cancellation of systematic uncertainties inherent in the use of (^*, we calculate 
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FIG. 2. Distribution of M^e in simulated events within the signal region overlaid with the 51 events 
observed in data. The interpretations of the various simulation histograms are identical to those 
in Fig. [TJ 



Ree by weighting each R\^ by the expected number of Dalitz decays, which is proportional 
to This gives 

We consider systematic uncertainties of the signal yields and efficiency ratios in our 
measurement of R^e- The systematic uncertainty in contributes a fractional uncertainty 
of 10.6%. Systematic uncertainties in yl^ and statistical uncertainties in ^* contribute a 
total fractional systematic uncertainty of 2.0%. Two systematic uncertainties remain; first, 
an uncertainty due to the different selection criteria on Mbc and 5M, as described above, 
which is estimated to be 4.1%; second, the uncertainty in the ratio of reconstructing an 
e^e~ pair to that of a 7. The uncertainty in this ratio is estimated by studying the decay 
ip{2S) J/ip 7r°7r°. The ratio between the number of events where one of the 7r° decays 
to 7e^e~ and the number of events where both vr" decay to 77 must be equal to twice the 
ratio of branching fractions -R^g = B{n'^ — )■ 7e+e~)/i3(7r° — )■ 77). Using this relationship, we 
measure Rf^ = (1.235 ±0.051)% in a manner similar to Ree by reconstructing the J/ ip vr^vr^ 
through these two decay modes of the 7r°. We restrict the energies of the to the range 20 to 
144 MeV (the mass difference rriD* — ''^dJ- Compared to the PDG value of (1.188±0.035)% 
[2I, we estimate the fractional systematic uncertainty on the ratio of efficiencies to be 6.5%. 
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FIG. 3. Distribution of Mbc of D*, 



"BC 

D^j events where Df — t- K'^ K^tt'^ . The points with 
error bars are data, the unshaded region is the D*'^ — )• Df^ signal. Background events in the 
highest shaded region (yehow onHne) are events with a direct paired with a photon from the 
D*~^ —7- Df^ decay. Background events in the middle shaded region (green online) are events with 
a direct Dj paired with a photon that did not come from -D*"*" — )• Dfj. Background events in the 
lowest shaded region (blue online) are from all other sources. 



These fractional systematic uncertainties are combined in quadrature to yield the final result 
for the ratio Ree'- 

Ree = [0.72+°:}^(stat) ± 0.10(syst)]%. (4) 

We use this result to re-evaluate the absolute branching fractions of the D*"*" meson as 
presented in Table IITI 

In summary, we report the first observation of a third decay mode of the -D*"*", the 
D*'^ — )■ D^e~^e~. We observe 51 candidate events in our signal region with an expected 
background of 12.6 events. The signal significance is 5.3a. The ratio of branching fractions 
Ree is measured as presented in Eq. and found to be consistent with our theoretical 
prediction. This implies that existing estimates of D*^ branching fractions should be revised. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, the Natural Sciences and 
Engineering Research Council of Canada, and the U.K. Science and Technology Facilities 
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TABLE II. Branching fractions in percent for the known decays of the D*^ meson from PDG 
2010 01, which assumed B{D*^ — D+e+e^) = 0, and our re-evaluation using the value of Ree 
reported in this Letter. 



Decay Mode 


PDG 2010 


This Analysis 


D*+ ^ D+j 


94.2±0.7 


93.5±0.5±0.5 


D*+ D + TT^ 


5.8±0.7 


5.8±0.4±0.5 


D*+ D+e+e- 





0.67t[J:}| ± 0.09 
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